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ABSTRACT 

We investigate the kinematics of ionized and neutral gas in a sample of 35 damped Lya systems 

flj . (protogalaxies) using accurate quasar spectra obtained with HIRES, the Echelle spectrograph on the 

^^ ' Keck I 10 m telescope. Velocity profiles with resolution of ~ 8 km s~^ are obtained for high ions such as 

C IV and Si IV, and for intermediate ions such as Al III. Combining these profiles with similar quality 

profiles obtained previously for low ions such as Fe II, we investigate the kinematic state of damped Lya 

protogalaxies in the redshift range 1.8 < z < 4.4 by comparisons between data for various ion pairs. 

We find the damped Lya protogalaxies comprise distinct kinematic subsystems: a low ion subsystem in 
which the low ions are physically associated with intermediate ions, and a high ion subsystem containing 
QQ I neither low nor intermediate ions. The evidence for two subsystems stems from (a) differences between 

f^ I the widths of the velocity profiles, (b) misalignment in velocity space between the narrow components 

0^ ' that make up the profiles in each subsystem, and (c) significant dissimilarities between the mean velocities 

^D , of the high ion and low ion velocity profiles. In every case we find that test statistics such as velocity 

^^ ■ width and various asymmetry parameters distribute differently for low ions than for high ions. We 

also find the absence of intermediate and low ions from the high-ion subsystem to indicate the latter is 
optically thin at the Lyman limit. 
^ i' Despite misalignment between their velocity components, the low and high ion kinematic subsystems 

are interrelated. This is indicated by detection of a statistically significant C IV versus low ion cross 
correlation function. It is also indicated by a systematic effect where the C IV velocity widths are greater 
than or equal to the low ion velocity widths in 29 out of 32 systems. These phenomena are consistent 
^ , with the location of the two subsystems in the same potential well. 

Subject headings: cosmology — galaxies: evolution — galaxies: quasars — absorption lines 

•i-H 

i^ , 1. INTRODUCTION component tends to occur at the profile edge. In PWl 

rrt T7.1, . 1 I. u j-j-j-i, and PW2 we used these properties to test various models 

. >T 1 ror the past several years we have been studying the r , , t , ttt i , i n^^ ; ^ i ; i 

kinematics of gas in damped Lya systems. We have fo- f°^" ^^"^P^^ ^ya systems. We adopted Monte Carlo tech- 

cused on the neutral gas (Prochaska & Wolfe 1997; 1998: ^^q^es by sendmg random sightlmes through gaseous con- 

u ti- Tnir"r~rTi77o\ 1 r- • j ;• ;i- hgurations spccifaed by the cosmological modci, the gcom- 

hereafter PWl, PW2) because of evidence suggestmg this '^ ^ ^ ■ c ^^ n ^- _l ttt ^i 

,1,1 ru c I. ■ 4.1- etry and physics of the configuration, etc. We then com- 

to be the source of baryons tor stars m current galaxies , i i- , -i . • c. ;..•.• i ,1 m i- 

t\xT ^c innc; i nn'7 t^ cc ic\c\a\ \xT j j-i, tittdt:^c putcd distributions ot test Statistics such as the prohle line 

(Wolfe 1995, 1997; Kauftmann 1996). We used the HIRES •1,1 1 1 ,1 -.i ,1 ..11., .1 

TT 1 1, , 1, nr i- inn/f\ -ti, t^ 1 T 1 n j. 1 Widths and compared them with the empirical distribu- 

Echelle spectrograph I Vogt 1994) on the Keck I 10 m tele- ,. tut n ... 1 • 1 1 i- i- 1 1 1 

,,,.,,•, ni r^ ■ i tious. We hrst tcstcd a Simple modcl lu which dark-mattcr 

scope to obtain accurate velocity pronles ot low ions such , , , , , . , . , , , . , 

fl^ Fe IT Si IT Ni IT and Al IT as thev trace the kmemalics ^^^^^ ^^1°^ enclose identical randomly oriented exponen- 

do PC 11, 01 11, IM 11, dllLL 1^1 11 db Liltjy LI dCt; Lilc KilicllldLlCo • l T 1 -.i T^ l ta r\r-r\ 1 1 1 

r , 1 ,1 rpi -1 r J-U' • -1-T, 1 Al tial disks With rotatiou speeds Vrot ~ 250 km s and 

of the neutral gas. ihe evidence tor this is the large Al . , , , ■ , . , „ on i n ■ i i- i 

H/Al HI ratios detected in most damped Lya systems, exponential scale-height h « 0.3i?,, where i?^ is the radial 

u-u-j-ij-i, • 1- -J • 4.U -ij scale-length (PWl). We found this model to satisfy all the 

which indicate the singly ionized species to be associated ..•.-,.. t. • • i i- i . ^ • •. 

• , 1 t- ^ ,-u t-i, ■ ■ A -A A i-u ■ ■ ■ statistical tests. Its principal disadvantage, however, is it 

with neutral rather than ionized gas, provided the ionizing . . , • , ttt i i • 

radiation is sunnlied hv external sources CProchaska and l*^ ^°* ^®* ^^ ^ COSmologlCal context. We then tested Semi- 

IdOIdLlOll lo oUUnilCQ Uy cXLt-llldi bUUiCco IX^iUClldoKd dllll i , • • r i i T i • /-^t^tv t ^ 

W If IQQfi'i analytic versions ol standard adiabatic CDM cosmologies 

mn , . 1 -i ci i-1 i- 1 J (SCDM) in which the neutral gas is confined to centrifu- 

ifie low-ion velocity profiles tfiat we fiave measured ^ ,, ' , i. , . , , ,, , , i r 

,, • , J- iu i 4. gaily supported disks m dark matter halos drawn from a 

comprise multiple narrow components tfiat are not ran- o ./ t-t- ,. ■ m, i • r i ■ 

,,,.,.,,,. 1 ■, Ti ii ii i i computed mass function, ihe predominance ot obiects 

domly distributed m velocity space. Rather the strongest .,, , ,, ,, i. ■ t • i.i .i . 

with low-mass and low Vrot, results m line widths that are 

^Visiting Astronomer, W.M. Keck Telescope. The Keck Obser- ^^^ Jq^ to match the observed distribution which extends 

vatory is a joint facility of the University of California and the Cali- . onr\ i —It 4' i j-i o/-<t^at i i 'j 

fornia Institute of Technology. *« 300 km s ^ . In particular, the SCDM models consid- 

ered by Kauffmann (1996) were found to be incompat- 
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ible with the kinematic data. Furthermore, Jedamzik & 
Prochaska (1998) showed that no CDM cosmology is con- 
sistent with models assuming the damped Lya systems to 
be single centrifugally supported disks within dark-matter 
halos. On the other hand Haehnelt et al. (1998) numer- 
ically simulated SCDM models with gas and dark mat- 
ter, and found that at z ^ 2 the damped Lya gas was 
distributed in low-mass protogalactic clumps rather than 
centrifugally supported disks. The combination of infall 
and chaotic motions apparently reproduce the kinematic 
data. 

This paper focuses on the kinematics of the ionized gas. 
Our aim is to place further constraints on galaxy formation 
models by studying gas located outside the neutral zones 
giving rise to damped Lya lines. Ionized gas is a generic 
feature of such models because initially the baryons are 
heated and ionized as they virializc in the potential wells 
of dark-matter protogalactic halos. However, the indicated 
velocity field is not unique, with some models predict- 
ing radial collapse to neutral disks (Mo & Miralda-Escude 
1996) , while others envision chaotic motions of ionized 
and neutral blobs (Haehnelt et al. 1997). As a result, de- 
termining the velocity structure of the ionized gas should 
help to clarify crucial events in the galaxy formation pro- 
cess (see accompanying paper by Wolfe & Prochaska 2000; 
hereafter paper II). 

In order to trace the velocity structure of the ionized 
gas we used HIRES to obtain accurate velocity profiles 
of the C IV, Si IV, and Al III ions; these are shown in 
§ 2. Throughout this paper we assume the C IV ion to 
represent the highly ionized gas, since we have a larger 
number of accurate C IV than Si IV velocity profiles. In 
§ 3 we construct frequency distributions of the profile line 
width, Av, for each ion. We compare these with each other 
and with the Af distribution for low ions such as Fe II or 
Si II which presumably trace the kinematics of the neutral 
gas. To further intercompare the low ion and high ion 
gas we examine correlations between the kinematics of low 
ions and high ions. To that end we consider the difference 
between the means of low ion and C IV profiles. We also 
consider the ratios of the Aw's for various ion pairs. We 
then cross-correlate the velocity profiles of various ionic 
pairs. In § 4 we present model- independent conclusions 
following from our results. 

2. DATA: OBSERVED VELOCITY PROFILES 

Table 1 lists the sample of damped Lya systems for 
which we have obtained velocity profiles. The spectra 
were acquired with HIRES at a resolution with FWHM 
« 6— 8 km s^^ and were reduced according to procedures 
outlined in PWl and PW2. The coordinate name of the 
background QSO is in column 1. Column 2 specifies the 
absorption redshift of the damped system. The entries in 
columns 3—5 contain the specific C IV, Si IV, and Al HI 
transitions used in our analysis. Entries with horizontal 
lines indicate that statistically significant profiles were not 
observed. Column 6 lists the low ion transitions used in 
our analysis and Column 7 gives the data references. Here 
and throughout the paper, the term low ion refers to ions 
such as Fe II and Si II, the dominant ionization states in 
neutral gas with large optical depths at the Lyman limit 
, i.e., tll > iC. The ionization potential (IP) of such 
ions is greater than 1 Ryd, i.e. 13.6 eV, while the IP of 



Table 1 
IONIC TRANSITIONS IN DAMPED Lyq SYSTEMS 





^abs 




Tr 
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QSO 


(J IV 


Si IV 


Aim 


Low 


Ref. 


QOlOO+1304 


2.309 


1548 


1393 


1854 


Ni II 1741 


W 


Q0149+33 


2.14075 


1548 


1393 


1854 


Fe II 1608 


W 


Q0201+3634 


2.4628 


1550 


1393 


1862 


Si II 1808 


W 


Q0216+0803b 


2.2930 


1550 


— 


1862 


Si II 1808 


S<^ 


Q0347-3800 


3.0247 


1548 


1393 


— 


Fe II 1608 


w 


Q0458-020 


2.03955 


— 


— 


1854 


Cr II 2056 


w 


Q0528-2505Ad 


2.14104 


1550 


— 


— 


Si II 1808 


s 


Q0841-0203A 


2.374518 


1548 


— 


1854 


Ni II 1741 


w 


Q0841-0203B 


2.476219 


1548 


1393 


1854 


Fe II 1608 


w 


Q0930+28 


3.23525 


1548 


1393 


— 


Fe II 1608 


s 


Q0951-04A 


3.85669 


— 


1393 


— 


Si II 1526 


w 


Q1055+46 


3.3172 


1548 


— 


— 


Fe II 1608 


s 


Q1104-18 


1.661375 


1550 


— 


1854 


Si II 1808 


s 


Q1108-07 


3.607619 


1548 


1402 


— 


Fe II 1608 


w 


Q1202-0725 


4.38285 


1548 


1393 


— 


Si II 1304 


s 


Q1215+3322 


1.9991 


1548 


1393 


1854 


Si II 1808 


w 


Q1223+17 


2.466083 


1550 


1402 


1862 


Si II 1808 


s 


Q1331+1704'' 


1.77636 


1550 


— 


1854 


Si II 1808 


w 


Q1425+6039 


2.8268 


1550 


1402 


— 


Fe II 1608 


s 


Q1759+7500 


2.6253 


1550 


— 


1854 


Si II 1808 


w 


Q1850+40 


1.99016 


— 


— 


1862 


Zn II 2026 


s 


Q1946+7658A 


1.7382 


1548 


— 


1854 


Si II 1808 


s 


Q1946+7658B 


2.8443 


1550 


1402 


— 


Si II 1304 


s 


Q2206-1958A'' 


1.920 


1548 


— 


1854 


Ni II 1741 


w 


Q2206-1958B 


2.07623 


1550 


1402 


— 


Al II 1670 


w 


Q2212-1626 


3.6617 


1548 


1402 


— 


Si II 1304 


s 


Q2230+02 


1.858536 


1550 


1402 


1862 


Si II 1808 


w 


Q2231-0015 


2.06615 


1548 


1393 


1854 


Si II 1808 


W-l-S 


Q2233+13 


3.14927 


1548 


1402 


— 


Fe II 1608 


s 


Q2237-0608 


4.0803 


1548 


1393 


— 


Al II 1670 


s 


Q2343+12 


2.42969 


1550 


1393 


1854 


Si II 1808 


s 


Q2344+12 


2.537789 


1548 


1402 


— 


Al II 1670 


s 


Q2348-1400'' 


2.2794 


1548 


1393 


1862 


Fe II 1608 


w 


Q2359-0203A=' 


2.095067 


1550 


— 


1854 


Si II 1808 


w 


Q2359-0203B 


2.153934 


1550 


— 


1862 


Si II 1526 


w 



"Data collected by our group 

''C IV profiles exhibit some saturation 

"^Data kindly provided by W. L. W. Sargent and collaborators 

''The 2 = 2.811 damped system toward Q0528— 2505 was omitted 

because its redshift exceeds that of the background QSO to which it 

may be associated. 



the next lower state (typically the neutral atom except for 
O) is less than 1 Ryd. The term high ion refers to ions 
where the next lower state has IP ^ 1 Ryd, for example 
Si IV and C IV. Photoionization studies show these ions 
to be in gas which is optically thin at the Lyman limit or 
in gas with moderate Lyman-limit optical depths (tll < 
10^). Finally, we introduce the term intermediate ions to 
describe those ions where the next lower ionization state 
has IP > 1 Ryd (e.g. Al HI). The same photoionization 
studies show that at column densities required for detec- 
tion (i.e. above lO^'^cm"^) these ions arise in gas which is 
optically thick at the Lyman limit. As a result, associa- 
tion of high ions with intermediate ions indicates the gas 
is optically thick at the Lyman limit. 

The velocity profiles are shown in Figures 1—5 which 
plot flux versus velocity, where the flux, I{v), is normal- 
ized to unit continuum. In Figures 1— 3 we compare C 
IV, Si IV, and Al III with the corresponding low ion pro- 
flies. The dark curves depict the high and intermediate 
ions and the grey curves depict the low ions. Figure 4 
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compares the C IV (dark curves) and Si IV (grey curves) 
profiles, and Figure 5 compares C IV (dark curves), Al 

III (light curves), and the low ions (grey curves). In all 
cases, V = km s^^ corresponds to the redshifts in Table 
1. The low ion transitions were selected on the basis of 
criteria outlined in PWl and PW2; i.e., for high signal- 
to-noise ratios and absence of saturation. The high ion 
transitions were selected according to the same criteria, 
where possible. In some cases the only transitions available 
exhibit saturation over sizable velocity intervals. The C 

IV profiles toward Q0216+08B, Q1331+17, Q2206-19A, 
Q2348— 14, and Q2359— 02A contain significant regions in 
velocity space with strong saturation. The same is true 
for the Si IV profile toward Q2348— 14. As emphasized 
in PWl, some of the statistical tests used to compare the 
data with model predictions are sensitive to saturation; in 
particular determination of test statistics measuring the 
profile asymmetries can be inaccurate in the presence of 
saturation. However, the results of statistical tests de- 
scribed below did not change qualitatively when the satu- 
rated C IV profiles were excluded from the sample. As a 
result the saturated profiles shall be included in all subse- 
quent analysis. 

The velocity profiles leave the following impressions: 

• In common with the low ion profiles, the intermediate 
and high ion profiles exhibit a multi-component structure 
comprising several narrow components spanning a wide 
range of velocity intervals. 

• The C IV and low ion profiles in Figure 1 are kinc- 
matically disjoined. In many cases strong C IV absorption 
components are at velocities where low ion absorption is 
weak or absent. In other cases C IV absorption is weak or 
absent at velocities where low ion absorption is strong. In 
many cases the velocity centroids of the C IV and low ion 
profiles appear to be significantly different. 

• Despite their differences, the C IV and low ion velocity 
profiles overlap in velocity space in such a way that the low 
ion profiles generally lie within the high ion profiles. 

• The low ion profiles also appear to be kincmatically 
disjoined from the Si IV profiles (see Figure 2). 

• The Al III and low ion profiles appear to be strongly 
correlated (Figure 3). 

• The C IV and Si IV profiles appear to be strongly 
correlated (Figure 4). 

• In many cases C IV absorption occurs in regions of 
velocity space that are free of Al III and low ion absorption 
(Figure 5). 

• In most cases the velocity widths of C IV and Si IV 
velocity profiles exceed the velocity widths of the Al III 
and low ion profiles. 

• Although the low ion and high ion velocity profiles 
exhibit "edge-leading" asymmetries in which the strongest 
velocity components are at the profile edges, the high ions 
tend to exhibit peaks at both edges more often than the 
low ions. 

In the sections that follow, we introduce several tests to 
quantify these phenomena. 



3.1. 



3. TEST STATISTICS 
Distributions for single ions 



To describe the kinematics of the gas in a quantitative 
manner, we consider test statistics that characterize the 



extent and shape of the velocity profiles. We discuss dis- 
tributions of such statistics drawn from the profiles cor- 
responding to given ions and then compare them. These 
empirical distributions are used to test theoretical models 
in paper II. 

In PWl we characterized the extent and asymmetry of 
the velocity profiles with four test statistics. Briefly stated 
they are: 

(1) Av, the velocity interval enclosing the central 90% 
of the integrated optical depth, jT(v)dv, where t{v) = 
\n[l/I{v)]. 

(2) fmm, the "mean-median statistic" given by \vmedian- 
Wmean |/(Aw/2), whcrc Vmedian and Vmean are thc median 
and mean velocities of the profiles. 

(3) fedg, the "edge statistic" given by |wpeafc-Wmea«|/(Aw/2), 
where fpeafe is the velocity of the absorption component 
with peak optical depth. 

(4) ftpk , the "two-peak statistic" given by 

±|wtpfc - f mean |/(Aw/2), whcrc Vtpk is the velocity of the 
component with second strongest peak optical-depth: ftpk 
is negative if Vpeak < Vtpk < Aw and positive if < Vtpk 
< Vpk- If necessary the optical depth profiles have been 
reflected in velocity space in such a way that fpeafe is always 

located at W < Vmean- 

Distributions for each test statistic are illustrated in Fig- 
ure 6 for thc low ion, Al III, Si IV, and C IV transitions. 
In PWl and PW2 we focused on the low ion distributions. 
We used them to argue in favor of models in which the 
low ions were confined to rapidly rotating disks (rotation 
speeds « 250 km s^^ ) and to rule out several other mod- 
els. We return to this topic in paper II. Here we check the 
null hypothesis that low ion and high ion kinematics stem 
from the same process. We do this by computing PkS: 
thc Kolmogorov-Smirnov probability that each high ion 
or intermediate-ion test statistic is drawn from thc same 
parent population as the corresponding low ion statistic. 
Figure 6 shows that the null hypothesis cannot be dis- 
missed, with two exceptions. First, the null hypothesis is 
highly unlikely when comparing low ion versus C IV distri- 
butions of Av since Pks(Au)= 0.002, and when comparing 
low ion versus Si IV distribution of Av since Pks{Av) = 
0.021: the high ions do not exhibit as much power at low 
Av as the low ions, and more power than the low ions at 
Av > 200 km s^^ . Indeed, the low probabilities suggest 
the low ions and high ions are in two distinct kinematic 
subsystems. Secondly, the null hypothesis is unlikely in 
the case of low ion versus C IV distributions of ftpk be- 
cause PKsiftpk) — 0.035; i.e., the C IV profiles display 
more double peaked profiles than the low ions. On the 
other hand Pvisiftpk) = 0.23 in the case of low ion versus 
Si IV profiles. Since compatibility of the Si IV and low ion 
distributions of ftpk cannot not be nfled out, this is the 
only case in which the C IV and Si IV kinematics differ. 
We suspect that statistics of small numbers may be affect- 
ing the Si IV distribution, a possibility we will test when 
larger samples become available. 

3.2. Correlations Between Macroscopic Kinematic 
Properties of Ion Pairs 

We now turn to correlations between kinematic proper- 
ties of gas in different ionization states. First we consider 
distributions of thc following new test statistic: 



IONIZED GAS IN DAMPED Lya PROTOGALAXIES 



(5) 5v — Wmean(ion o) — Wmean(ion 6), i.e., the difference 
between the means of the velocity profiles of ion a and ion 
b. 

Figure 7a shows the distribution of 5v when a and b are 
(i) low ion and C IV, (ii) C IV and Al III, (iii) C IV and 
Si IV, and (iv) low ion and Al III. Not surprisingly, the 
distribution widths arc narrower f or C IV versus Si IV 
and low ion versus Al III than for low ion versus C IV. 
This just indicates that the high ion subsystem includes 
the high ions C IV and Si IV, but not the intermediate 
Al III ions which are mainly associated with the low ion 
subsystem. The standard deviations for the 3 distributions 
are given in Table 2. Although the dispersion of the low 
ion versus C IV distribution is relatively high, as — 67 



km 



it is sufficiently low to place crucial restrictions 



on most models (see paper II). 

Next we consider correlations between the Aw of vari- 
ous ion pairs. In Figure 8a we plot AwciV) the Aw inferred 
from the C IV profiles versus Awiow, the Aw inferred from 
the low ion profiles. In Figures 8b and 8c we do the same 
for Awciv versus AwsiiVj the velocity width of Si IV, and 
Awaiiii, the velocity width of Al III versus Awiow. The fig- 
ure shows Awciv to be correlated with Awsiiv and Awaiiii 
with Awiow. This is demonstrated in Table 3 showing the 
Kendall r correlation coefficients. The only cases with sta- 
tistically significant r (i.e., at the ha level) arc for the C 
IV versus Si IV correlation and the Al III versus low ion 
correlation. We conclude that statistically significant cor- 
relations exist between Awciv and AwsiiVj and between 
AwAiiii and Awiow. Evidence for such correlations within 
the remaining ion pairs is marginal. On the other hand 
Figure 8a shows that for all but 3 objects, Awciv > Awiow, 
and even the three exceptions nearly satisfy this relation. 
In other words the low ion velocity width acts as a "fioor" 
to the high ion velocity widths. While not a correlation, 
this is an important systematic effect indicating that the 
high ion and low ion systems are interrelated in a way that 
must be accounted for by models of the ionized gas. 

The relationship between the velocity widths is quanti- 
tatively described by the test statistic: 

(6) f ratio = Awa/Awfe, the ratio of the velocity widths 
of ion a to those of ion b. 

The corresponding distributions arc shown in Figure 7b. 
Notice how the f ratio for low versus C IV and C IV versus 
Al III are uniformly distributed at 0.7 < f Ratio < 3, while 
the C IV versus Si IV and low versus Al III distributions 
exhibit significant peaks at f ratio = 1. This is further 
evidence for an association between Al III and low ions, 
and between C IV and Si IV. 

We searched for evidence that other test statistics were 
correlated between the different ion pairs. None was found, 
with the possible exception of correlations between the Al 
III and low ion f^dg which exhibit Kendall t coefficients at 
the 3.0(7 level. Correlations between between the C IV and 
Si IV ftpk were found at the 2.5cr level. We also checked for 
correlations between pairs of test statistics corresponding 
to a given ion. Ledoux et al. (1998) reported evidence for 
a low ion correlation between f^dg and Aw out to Aw = 
150 km s~^ , but found no such correlation at Aw > 150 
km s"""^ . They claimed the "break" at 150 km s~^ argued 
against the presence of damped Lya systems comprised 
of disks with rotation speeds greater than 150 km s^^ . 



It is difficult to assess the validity of their claim since no 
quantitative estimate was given for the significance of the 
correlation. With our larger data set (39 versus their 26 
systems) we find evidence for such a correlation at the 3.6(t 
level, but no convincing evidence for a "break" at Aw « 
150 km s^^ . In fact when systems with Aw > 150 km s^^ 
were eliminated, the significance level of the correlation 
decreased to 2.8(7. This argues against the presence of such 
a "break" and is consistent with a correlation extending 
to Aw « 250 km s~^ . As a result disks with rotation 
speeds exceeding 150 km s^^ cannot be ruled out. In any 
case, the correlation between /e^g and Aw is tentative (see 
Figure 9), and needs to be tested with larger data sets. 

3.3. Cross Correlation Functions 

We next investigate whether the detailed component 
structures exhibited by the velocity profiles of various ions 
are correlated. In Figure 7c we plot ^ab{v), the cross- 
correlation functions between ions a and b versus lag ve- 
locity w, and the la errors, a^ ^(w). The ^ab{v) are given 

by 



Cah(w) 



1 ^ 

-T 



"fc 



fe=i ^/Km^) -1 



Tak{ui)][TbkiUi - w)] 



where the first sum is over N profiles, and the second is 
over Uk velocity channels, Ui, spanning the /c*'' absorption 
profile. The quantity Tak(v) is the optical-depth profile 
of ion a for the k*^ profile. The quantity, £,aai^)t is the 
autocorrelation function for ion a in the k*^ profile at w = 
0, and is defined as 



.k 

Saa 



iO)=Y.[rakiu^W 

i=l 



(2) 



As a result the cross-correlation functions are normalized 
so that perfectly correlated ions yield ^ab{v) = 1 at zero 
lag velocity. The errors are dominated by cosmic vari- 
ance and are computed by a "bootstrap" method (Efron 
& Tibshirani 1986). 

Comparison between the cross-correlation functions and 
the errors demonstrates that in some cases the measured 
£,abiv) are statistically significant out to w ^^ 150 km s^^ . 
The strongest correlations shown are between the low ion 
and Al III profiles and between the C IV and Si IV profiles. 
This confirms the rather accurate one-to-one alignment be- 
tween the velocity components comprising these ion pairs 
(see Figures 3 and 4). The 20 km s~^ half- width of ^ab 
for Al III versus low roughly equals the half widths of the 
wider components in Figure 3, which supports the idea of 
common velocity components for Al III and low ions. The 
30 km s~^ half width of ^ab for C IV versus Si IV argues 
for common component structures for these two ions. 

By contrast the £,abiv) for the C IV versus low ions or C 
IV versus Al III exhibit significantly lower, but still statis- 
tically significant, values at w = (the differences exceed 
6a). The lower ^ab(O) stems from misalignment between 
the low ion (or Al III) velocity components with the C 
IV components in many of the profiles. In this case the 
correlation amplitude arises from overlap in velocity space 
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between the ion pairs rather than one-to-one ahgnnient be- 
tween individual components. This interpretation is sup- 
ported by the ^ 80 km s~^ half widths of £^ab for the C IV 
versus Al III or C IV versus low profiles. These half-widths 
more closely resemble the coherence scales of contiguous 
multiple components in the C IV profiles than the half 
widths of individual components (see Figure 1). 

4. SUMMARY AND CONCLUDING REMARKS 

Using high-resolution absorption spectra obtained with 
HIRES, the Echelle spectrograph on the Keck I telescope 
(Vogt 1994), we have probed the kinematics of ionized gas 
in a sample of 35 high-redshift damped Lya protogalaxies. 
Specifically, we obtained velocity profiles of the high ions C 
IV and Si IV and the intermediate ion Al III. We studied 
the kinematic state of the gas by constructing empirical 
test statistics which characterize (a) the widths and sym- 
metry properties of high ion, intermediate-ion, and pre- 
viously obtained low ion velocity profiles, and (b) corre- 
lations between the kinematic properties of various ionic 
pairs. 

What have we learned that is new? We answer this ques- 
tion by discussing model-independent conclusions inferred 
from the data: 

(1) The damped Lya absorbers comprise distinct kine- 
matic subsystems: a low ion subsystem in which low ions 
such as Al II are physically associated with intermediate 
ions such as Al III , and a high ion subsystem consisting 
of ions such as C IV and Si IV. This is indicated by the 
similarity between the C IV and Si IV velocity profiles and 
between the low ion and Al III velocity profiles, and by the 
misalignment of the velocity components comprising the C 
IV and low ion profiles. 

(2) The low ion and high ion kinematic subsystems are 
related despite misalignment of their velocity components. 
This follows from the detection of a statistically signifi- 
cant C IV versus low ion cross-correlation function which 
exhibits lower amplitude and a wider half-width than ei- 
ther the low ion versus Al III, or C IV versus Si IV cross 
correlation functions. Whereas the high cross-correlation 
amplitudes of the latter two ionic pairs arise from the 
coincidence between corresponding velocity components, 
the lower amplitude of the C IV versus low ion cross- 
correlation function is due to a general overlap in velocity 
space between the line profiles. In any case the relation 
between low ion and high ion subsystems is further indi- 
cated by a systematic effect in which Auciv ^ Auiow in 
29 out of 32 profiles. 

(3) The difference between the mean velocities of the 
C IV and Si IV velocity profiles or between the low ion 
and Al III profiles exhibits distributions with dispersions 
equaling 23 and 30 km s~^ respectively. By contrast the 
low ion versus C IV distribution has a significantly wider 
dispersion of 67 km s^^ . This is more evidence for distinct 
kinematic subsystems. 

(4) The absence of intermediate ions and low ions from 
the high ion subsystem indicates the latter is optically thin 
at the Lyman limit. The lack of mixed ionization states 
distinguishes the high ion gas in damped Lya protogalax- 
ies from high ion gas in the ISM (Savage et al. 1993), high- 
redshift Lyman-limit systems (Prochaska 1999; Prochaska 
and Wolfe 2000), and z ~ 0.7 Mg II selected absorbers 
(Churchill et al. 1999) where the velocity profiles of the 



low ions and intermediate-ions resemble those of the high 
ions. 

These results have rather general implications. First, 
kinematic subsystems placed in the same potential well 
generate velocity profiles that tend to overlap in velocity 
space. In cases where the velocity field of the neutral gas 
is constrained to fewer degrees of freedom than the ionized 
gas the resulting profiles will be characterized by Auciv ^ 
Auiow. Therefore, we shall test the hypothesis that both 
kinematic subsystems are subjected to gravitational forces 
arising from the same mass distribution in paper II. Sec- 
ond, the absence of ionized gas optically thick at the Ly- 
man limit implies a distinction between the damped Lya 
systems on the one hand and the Lyman limit systems 
and Mg II absorbers on the other. This is potentially sig- 
nificant because most models assume that both classes of 
absorber arise in the same physical object. Finally, in pa- 
per II we use the empirical distributions of test statistics 
as well as the empirical C IV versus low cross-correlation 
function to constrain semi-analytic models of galaxy for- 
mation. Specifically we focus on models in which ionized 
gas in virialized dark-matter halos falls onto centrifugally 
supported neutral hydrogen disks. 

We wish to thank Eric Gawiser and David Tytler for 
valuable discussions. We also thank Wal Sargent for gen- 
erously giving us data prior to publication. AMW was 
partially supported by NSF grant AST0071257 and JXP 
acknowleges support from a Carnegie postdoctoral fellow- 
ship. 
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Fig. 1. — Velocity profiles of C IV (dark lines) and low ion (grey lines) transitions for all damped Lyo protogalaxies in Tablel in which 
HIRES C IV and low ion profiles have been measured. The velocity v = corresponds the absorption redshift, z^bsi shown in the Table. 
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Fig. 2. — Velocity profiles of Si IV (dark lines) and low ion (grey lines) transitions for all damped Lya protogalaxies in Table 1 in which 
HIRES Si IV and low ion profiles have been measured. Velocity scale same as in Figure 1 



Table 2 
DISPERSIONS OF 5v DISTRIBUTIONS 



Dispersion 

cr5t,'*(kni s^^ 



C IV versus Low C IV versus Si IV Al III versus Low 



67.0±8.2 



23.2±3.4 



30.0±4.2 cr 



"One sigma error in agy computed by assuming (n— l)s'^/o''^ is distributed as x^(n— 1) where n is the number of degrees of freedom, x^(n— 1) is 
X^ for n— 1 degrees of freedom, and where s^ and a^ arc the measured and true variance of ag^ (Frodesen et al. 1979) 



Table 3 
CORRELATIONS BETWEEN IONIC VELOCITY INTERVALS 

Kendall t 

Aim 



ION C IV 



"snv 



Low 



CIV 
Si IV 

Aim 



0.80±0.16 0.57±0.16 0.33±0.16 

0.76±0.22 0.52±0.15 

0.73±0.16 
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Fig. 3. — Velocity profiles of Al III (dark lines) and low ion (grey lines) transitions for all damped Lya protogalaxies in Table 1 in which 
HIRES Al III and low ion proflics have been measured. Velocity scale same as in Figure 1 
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Fig. 4. — Velocity profiles of C IV (dark lines) and Si IV (grey lines) transitions for all damped Lyo protogalaxies in Table 1 in which 
HIRES Si IV and C IV profiles have been measured. Velocity scale same as in Figure 1. 
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Fig. 5. — Velocity profiles of C IV (dark lines), Al III (light lines), and low ions (grey lines) transitions for all damped Lya protogalaxies 
in Table 1 in which Al III, C IV, and low ion profiles have been measured. Velocity scale same as in Figure 1 
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Fig. 9. — Plot of /eda versus Av for low ion profiles in Table 1 



